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In this paper we analyze the compatibility of two- zero neutrino Majorana textures with the recent 
experimental data. Differently from previous works, we use the experimental data to fix the values 
. of the non-vanishing mass matrix entries and study in detail the correlations and degree of fine- 

| tuning among them, which is also a measure of how naturally a given texture is able to describe 

all neutrino data. This information is then used to expand the textures in powers of the Cabibbo 
angle; extracting random 0(1) coefficients, we show that only in few cases such textures reproduce 
the mixing parameters in their 3<j ranges. 
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I. INTRODUCTION 

> 

reactor angle #13, as derived from the analysis of — > v e appearance and v e — > v e disappearance flavour transitions. 
In particular, Daya Bay and RENO provide a clear evidence of several standard deviations from zero of the reactor 
angle: 



The recent T2K Q, MINOS fl, Daya Bay Q and RENO Q| results contain a strong indication in favor of a large 



^1- 

^ . sin 2 20i3 = 0.092 ± 0.016 ± 0.005 (Daya Bay) sin^ 20 X3 = 0.103 ± 0.013 ± 0.011 (RENO) . (1) 



This has recently pushed the model building industry to study the possibility to accommodate a large 0i 3 , together 
with the well established solar and atmospheric mixing angles In spite of the good accuracy in the mixing 
parameter determination, the structure of the neutrino mass matrix is still unknown; two-zero textures of a Majorana 
type have been largely studied before the T2K and MINOS data [6] and an update analysis (including the hints for 
a large 6*13) has been recently presented in Q and Q. The common approach of such studies is that of using the 
zeros of the Majorana mass matrix to get simple relations among the neutrino mixing angles and mass differences, 
independently on the parameters defining the mass matrices. Such relations can then be used to determine the three 
neutrino masses and three CP-violating phases. It turns out that some of the two-zero textures are still compatible 
with the recent data at the 3tr level whereas all the textures marginally compatible or incompatible with the pre-T2K 
and MINOS results are now definitely ruled out 0]. In addition, a look at Fig. 11 of the same Ref.0] shows that also 
the texture called B4 (see later for its explicit form) has a problem with $13 , since its distribution at the 3a level does 
not touch the value 0i3 = 7° and is then marginally compatible with eq.([T]). 
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Although the procedure outlined above is very powerful to understand how the mixing angles and phases are related 
to each others (for a given texture), it gives no hint on the order of magnitude of the parameters defining the mass 
matrices; it also does not reveal whether the matrix elements are correlated and/or if they have to be fine-tuned 
(and to what degree) to reproduce the experimental data. In this letter we want to clarify all these issues and study 
in details the features of the mass matrix elements needed to accommodate the experimental results within a given 
texture. In principle one can use different approaches; here we rely on two different estimators. We first perform a 
X 2 analysis of the still allowed two-zero textures, with the aim of getting a global information on the goodness of the 
mass matrix to fit the neutrino data. For every texture, we explicitly show the values of the parameters that fit the 
data, in such a way that the reader can easily notice the deviation of the mass matrix parameters from the natural 
(in absence of a more fundamental flavour model) 0(1) value; such a deviation gives an idea of the naturalness of the 
textures. 

Beside the minimum of the x 2 ; the fit procedure will also return the corresponding values of 5 neutrino observables, 
namely the two independent mass differences Am 2 ol and Am^ im and the three mixing angles #13, #12 and #23- We 
do not include the Dirac CP phase into the x 2 definition since its value has not been measured yet. The minimum of 
the x 2 cannot be the only criterion to classify mass models, since a relatively good value can be obtained at the prize 
of a strong fine-tuning among the mass parameters. For this reason, we also give for each minimum a quantitative 
measurement (our second estimator) of the amount of fine-tuning needed. This is obtained calculating, for every 
parameter, the shift from the best fit value that changes the x 2 by one unit, with all other parameters fixed at their 
best fit values. 

Since our analysis is intended to make a relative comparison for the degree of fine-tuning between various textures, 
we do not take into account the larger fine-tuning implied by the tiny errors on the masses in the charged lepton 
sector, which would otherwise make our comparison completely useless. 

The results of such a procedure for a normal order (NH) of the neutrino mass eigenstates is described in Section [TT1 
Using the best fit values of the mass matrix parameters, we expand the textures in terms of the Cabibbo angle A = Ac, 
with generic complex coefficient with absolute values of 0(1) for any non- vanishing matrix elements. In this way, we 
automatically take into account the hierarchy among different entries of the mass matrix but it is not obvious that 
these coefficients combine to reproduce the best fit points of the neutrino observables. It is then interesting to study 
separately the degree of predictivity of such textures, extracting randomly the complex coefficients and analyzing the 
values of the obtained mass differences and mixing angles. This procedure and its outcomes are detailed in Section 11111 
while in Section llVl we repeat the same exercise for the inverted hierarchy. Section [V] is devoted to our conclusions. 



II. OBSERVABLES AND ANALYSIS DETAILS 



Although there exist a large class of Majorana mass matrices with two-zero entries, we restrict our analysis to the 
seven models still compatible with the data, according to Q, from which we also adopt here the same nomenclature: 
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A general complex symmetric matrix contains 6 moduli and 6 phases, for a total of 12 independent parameters. The 
redefinition of the neutrino fields allows to eliminate three phases; other two can be thrown away if they are associated 
with the vanishing elements of the mass matrix. In total, for two-zero textures we are left with 4 moduli and one 
phase ip that, as it can be easily demonstrated with a trivial computation, can always be associated to 04, for a total 
of 5 independent parameters. These are the quantities to be determined by the \ 2 fit- 

It has been pointed out in 7] that the zeros in A\ — A2 as well as those in B\ — B2 and B3 — B4 are related by a 
P23 symmetry 




Paa = 1 (5) 



in such a way that: 



and so on. This in turn implies: 



A 2 = P23 A x P 2 J 3 (6) 
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and similarly for the pairs B12 and B3 4. Analytical as well as numerical considerations on the range of the mixing 
parameters implied by the previous textures have been already analyzed in the literature (see for example Q); the 
main conclusions are that the textures Ai describe the data much better in normal hierarchy than in the inverted one 
whereas the texture C allows for the inverted hierarchy. For the B^s, the mass spectrum is quasi-degenerate and it 
depends on the assumed octant of #23; for #23 < 7r /4 -81,3 are compatible with the normal hierarchy and 52,4 with 
the inverted one (the behavior is opposite for 623 in the second octant). 

In order to quantify the qualitative behaviors described above, we start with the normal hierarchy (NH) case, using 
the best-fit points and the 1<t uncertainties on the mixing parameters from 9], summarized in TabJH The definition 
of the x 2 is as follows: 



[sin 2 ^-(sin 2 %)^] z | N-(A^H^ 

X —^^' = 12,13,23 2 <~ i -'k=sol,atm § ' \°) 



where (sin 2 9ij) data and (Am|) da * a are the values of the mixing angles and mass differences reported in Tab|H re- 
spectively, and a the corresponding ler errors also quoted in Tab|l] sin 2 0^ and Ato 2 depend on the five unknowns 
ai ; .. ; 4 and tp. Given that two statistically equivalent best fit values for #23 emerged from the analysis in Q, for every 
textures we fit both values and retained the ai,..,4 and ip corresponding to the minimum among the two % 2 . 

To estimate the degree of fine-tuning we used the parameter cLft defined in [12j. This adimensional quantity is 
obtained as the sum of the absolute values of the ratios between each parameter and its " error" , defined for this 
purpose as the shift from the best fit value that changes the % 2 by one unit, with all other parameters fixed at their 
best fit values (this is not the error given by the fitting procedure because in that case all the parameters are varied 
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0.59 0.023 


la 
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Best fit (IH) 


7.50 
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Iff 
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[-2.39, -2.50] 
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TABLE I: The latest global-fit results from 
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at the same time and the correlations are taken into account): 



Aft = ^ 



pan 



err. 



(9) 



It is clear that dpx gives a rough idea of the amount of fine-tuning involved in the fit because if some \erri/pari\ 
are very small it means that it takes a minimal variation of the corresponding parameter to make a large difference 
on the x 2 ■ As in [l2| , we can compare this quantity with the sum of the absolute values of the ratios between each 
observable and its error as derived from the data: 
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that for the set of data in Tab|T]is doat 
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We have summarized our numerical results in TabHTl where we have shown the value of the Xmim the fine-tuning 
parameter dpp, the best fit values of the neutrino matrix elements and the predictions of some of the corresponding 
experimental observables. Notice that we did not display the values at the best fit (cii,ip) for 612 and the mass 
differences, since we always got values indistinguishable from Tab|TJ First of all, we see from Tablrfl that for all 
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TABLE II: Fit results. Shown are the best fit values of the five independent parameters defining the two-zero textures as well 
as the obtained values of the mixing angles, the \ 2 an d the fine-tuning dFT of the minimum. Results refer to the NH case. 

patterns but C we can find a very small \ 2 ! this is not particularly surprising since, as claimed in the Introduction, 
all patterns are compatible with the experimental data and C prefers the inverted hierarchy. Differences arise as we 
deal with the chosen estimators, the magnitude and hierarchies of the matrix elements and the values of the mixing 
parameters corresponding to the best fit points. In fact, patterns Ax^ show good agreement of all observables with 
Tab H] and the fitted parameters (matrix elements) are close to 1, with only modest dpT values. On the contrary, 
although the textures 5^=1,. ..,4 also have small x 2 's, the fine-tuning among the matrix elements is quite strong, as 
it can be seen from the larger dpT values compared with the A\^ ones. In addition, a more pronounced hierarchy 
between the parameters is needed to fit the data. The Bi textures only differ for the preferred #23; the reason can be 
traced back to eq.©: B13 prefer an atmospheric angle less than 7r/4 and very close to the best fit sin 2 6*23 ~ 0.41 
whereas textures i?2,4 would rather favor a 6*23 > 7r/4, very close to the other best fit sin 2 6*23 ~ 0.59. Notice that 
these results are slightly dependent on the set of data used for the fit; if, instead of TablH we had adopted the best 
fits and errors as given in (lCf . with only one minimum at sin 2 $23 ~ 0.52, the situation would have been different 
with, for example, a better \ 2 an d smaller dpr for B2 compared to B\. This would also be the case for A\^ but, 
since these textures do not prefer any octant, both matrices perform very well. The case C is the more finetuned and 
it prefers #23 = 7r/4, that is partially excluded with the current data. 

Summarizing, textures Bi and C can be classified as less natural with respect to A\^ because of the larger fine- 
tuning parameter dpT and the more pronounced hierarchies among the <Zj matrix elements. It is important to stress 
that, given the relatively large number of fitted parameters, for some textures several good x 2 comparable in size 



■5 



can be found, with matrix elements slightly different from fit to fit. In these cases, we observed that the smaller % 2 
is usually given by very fine-tuned solutions with a strong hierarchy between the and, as a selecting criteria, we 
decided to present the results with the smaller value of the product % 2 • dpx- 



III. TEXTURE PREDICTIVITY 

From the model building point of view, it could be interesting to describe the previous textures in terms of powers 
of the Cabibbo angle A. Taking the values of the ai coefficients in Tab|n]and maintaining the hierarchies among them 
we get: 
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B 4 : X 2 d 2 Xdi a 3 , C: Aa 2 a 4 , (13) 



where the new parameters ai are unspecified 0(1) numbers. In the spirit of U(l) constructions, the coefficients of 
every matrix elements are unspecified complex entries with absolute values of 0(1); then, we have not expressly shown 
any possible correlation among the matrix elements, although they can be important to get the correct values of the 
neutrino mixing parameters. The only thing we can say is that there exists a non-vanishing volume in the parameter 
space where the neutrino data can be fitted with sufficient accuracy. 

Some analytical considerations on the predicted masses and mixing angles are possible by means of standard 
perturbation theory in the expansion parameter A applied to eqs. (|imi3| . It turns out that A\ >2 predict almost 
maximal # 23 and a value of the reactor angle as large as A; on the other hand, at leading order (LO) 9\ 2 is vanishing 



and receives corrections at 0(A). Textures of structure similar to -Bi, 2j3 ,4 have been carefully studied in [ll( where, 
however, some of their many properties have been obtained in the hypothesis that the matrix elements have well- 
definite values and not unspecified 0(1) coefficients in front of them. For vanishing A, the -Bi j2 , 3j4 have two degenerate 
eigenvalues; if the third one is the largest, then we can associate the two degenerate masses to the solar sector and 
(1,0, 0) T to m 3 . When also the corrections from A are taken into account, the previous assignments produce small 
r = Am 2 oi /Am 2 tm and very large #i 3 . The solar and atmospheric angles are unstable and can get almost any value. 
On the other hand, if the non-degenerate mass is smaller than the degenerate pair of eigenvalues, then we are forced 
to assign such a smaller mass to mi and the other two eigenvalues to m 2 and rn 3 . Again, reintroducing A, we get 
large r, small 0i 2 ,i 3 ~ 0(A 2 ) and the correct value 6* 23 = 7r/4 + 0(A). This last case seems to agree better with 
the experimental results (as it can be seen from the values of the at in Tab|TI]) although a moderate finetuning on 
the parameters is needed (see later for the correlations). Texture C has the same LO as for B so, again, we have to 
distinguish two cases as before. When the almost degenerate eigenvalues are associated to mi,2, then r is vanishingly 
small, $i 3 , 23 ~ 0(1) an d the solar angle is always small, suppressed by the (11) entry of the neutrino mass matrix. 
In the opposite case, r ~ 1, #i 3i i 2 ~ 0(A) and 6* 23 ~ 7r/4. Again this last case is favoured by the data (always at the 
price of a large finetuned correlations between the a,i). 



These considerations are confirmed by a numerical analysis; we performed a MonteCarlo simulation extracting 
randomly all 0(1) entries of eas. (1111115)) with moduli in [1/3, 3] (with equal probability of obtaining a number smaller 
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than 1 and larger than 1); the unique phase (p is extracted flat in [— 7r,7r] and can always be associated to 04 (ljj ; 
we then built the hermitian matrix mj,m y and extracted the eigenvalues and eigenvectors. We do not impose any 
external additional constraints on the mass differences and mixing angles. Notice that this is different from what 
studied, for example, in p| since we started from explicit mass textures (derived from our fit procedure) and use a 
typical f (1) approach to estimate the predicted angles whereas the authors of [7| generate sets of random numbers of 
experimentally allowed (#12, #13, #23) an d mass differences to predict other quantities of interest and retaining those 
mixing parameters satisfying certain consistency relations. We plot the distributions of the ratio r = A?7j 2 o; /Am 2 (m 
and sin 2 0i3, sin 2 023 and sin 2 812 in Figs J1I2I3I and l4l respectively. Given the 23-symmetry discussed around eq.Q, 
which relates the mixing angles among those textures connected by P23, we do not show the obtained distributions 
for all textures but rather prefer to superimpose in the left panel the distributions for the textures B\ and B3 (solid 
and dashed lines, respectively) and in the right panel the results for A\ and C (dashed and solid lines, respectively). 
The vertical black lines indicate the 3cr allowed region on the variable under discussion. 





FIG. 1: Distribution for the variable r as obtained from the textures B\ and B3 (left panel - solid and dashed lines, respectively) 
and Ai and C (right panel - dashed and solid lines, respectively). The vertical black lines indicate the 3a allowed region on r. 
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FIG. 2: Distribution for the variable sin 2 $13 as obtained from the textures Bi and B3 (left panel - solid and dashed lines, 
respectively) and Ai and C (right panel - dashed and solid lines, respectively). The vertical black lines indicate the 3a allowed 
region on #13. 



As expected, it is difficult to reproduce the small value of r ~ 0.03 in a natural way, since it usually requires a 
large amount of fine-tuning on the mass parameters. We observed only a moderate depletion of r for the patter A\, 
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FIG. 3: Distribution for the variable sin 2 823 as obtained from the textures Bi and B3 (left panel - solid and dashed lines, 
respectively) and Ax and C (right panel - dashed and solid lines, respectively). The vertical black lines indicate the 3a allowed 
region on 823- 




sin 2 6>i2 sin 2 #i2 

FIG. 4: Distribution for the variable sin 2 thetax2 as obtained from the textures Bx and B3 (left panel - solid and dashed lines, 
respectively) and Ax and C (right panel - dashed and solid lines, respectively) . The vertical black lines indicate the 3a allowed 
region on 812 . 

r ~ 0.5 (FigO}. From Fig(2]we see that A\ is very powerful in reproducing the relatively large values of #13, a feature 
not shared by the textures Bx t 3 since they present a peak in the distribution at the lower bound of the 3<7 range. 
Textures C performs worse than the others, since even the smallest peak at sin 2 6* 13 ^0.1 is essentially excluded. We 
take this as an indication that the texture C is not particularly suitable to fit the data in normal hierarchy. The 
angle #23 (Fig{3]) is well reproduced by essentially all textures, with deviations as discussed below eq. (fl"3"|) . Also, the 
distribution of 612 is almost flat for and favors small values for A\ and C (Fig|4|. 

It is interesting to evaluate the fraction of points which simultaneously give 9%j and r into their the 3tr bounds. 
This can be considered as another estimator of the fine-tuning among the matrix elements needed to reproduce the 
correct experimental data. The computation results in the following small fractions: 

{Ax, A 2 , B x , B 2 , B 3 , B^ C} = {0.1, 0.1, 0.004, 0.003, 0.003, 0.002, < 0.0001}% , (14) 

confirming that the textures Ax,2, derived from Ax, 2 of eq.©, are the most appropriate to describe the whole data in 
the neutrino sector in NH. 

The previous plots do not give any indications of the degree of correlations among the matrix elements of eas. dllt - 
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[T3f needed to get the desired values of the mixings and r. This is illustrated in Figj5l where we show, for all seven 
textures, the correlations among ai and 0,2,3, |&4 1 - 
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FIG. 5: Correlations among 2i anci 52,3, | S4 1 /or t/ie textures studied in the paper, for the case of NH. From above to below: 
Ai, A2, Bi, B2, -B3, B4, C . 
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The results of Figj5]show that the textures Ai reproduce the data at the price of mild correlations between all the 
parameters, whereas we always found a significant correlation between a pair of hi (the ones corresponding to the 
largest in Tab Hi]) for B^ The strong correlation observed for the C is a clear signal that such a texture is specially 
suitable to describe the data in IH more than in NH, so the matrix elements have to conspire in order to reproduce 
the data in the "wrong" hierarchy. 

We also observe that the degree of fine tuning for the Ai , Bi and C textures is the same as for the original matrices 
because numerator and denominator in the definition of cLft scale with the same power of A. 

IV. RESULTS FOR THE INVERTED HIERARCHY 

We have repeated the same fit procedure and analysis for the inverted hierarchy (IH). The main results are: 

• the textures A\^ have now a very worse Xmin ~ 0(100) and then are not suitable to describe the neutrino data; 
we do not discuss such cases any further; 

• the agreement of the texture C is improved with respect to the normal hierarchy case, both in terms of fine-tuning 
and correlations; the corresponding expression expanded in terms of A, C IH is given by: 

( l 1 *\ 

C IH : 10 1; (15) 
\1 1 0/ 

• the main features of the textures Bi remain essentially the same, with large correlations among the parameters, 
so that the expanded textures are the same as for the normal hierarchy case, that is: B{^ 3 4 = -61,2,3.4; 

For the sake of completeness we report such results in Tab llHI 
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TABLE III: Fit results for the inverted hierarchy. Shown are the best fit values of the five independent parameters defining the 
two-zero textures as well as the obtained values of the mixing angles, the \ 2 and the fine-tuning dpx of the minimum. 

The analytical discussion for the textures B IH can proceeds along the same lines as for the normal hierarchy case: 
in fact, the smallest mass m.3 can be associated to the eigenvector (1, 0, 0) T (thus giving small r and close to maximal 
013, with undetermined #12,23) or to (0,1,-1) (producing r ~ 0(1), vanishing #12,13 and maximal #23)- For the 
texture C IH we can get an analytical understanding of the predicted masses and mixing angles working with real 
coefficients close to 1; in this case, the texture predicts (in absence of accidental cancellations) large mixings for every 
angles #13,12,23 and quite a large mass ratio r. Except for the atmospheric angle, whose values are broadened in the 
whole [0, 7r/2] interval, these features are maintained as we allow for generic complex coefficients with absolute values 
in the [1/3,3]. 

In Fig[6] and FigJT] we have shown the distributions resulting from the numerical analysis in which the matrix 
elements of the tilded textures are extracted randomly with moduli in [1/3,3] and phases in [— 71", 7r]. Results are 



10 



shown only for B[ H (solid line), B^ H (dashed line) and C IH (thick dashed line). As it is common for inverted 
hierarchy, it is more difficult to reproduce the correct value of r since this needs a stronger conspiracy among the 
matrix elements. This is the reason why the distributions shown in the left panel of Fig(B] are more shifted toward 
r ~ 1 than in the NH case but for C IH , which gives very large 6*13. For the distributions of the other mixing angles, 
we do not observe any significant difference with respect to the NH case. 

Summarizing, there is no preferred patterns in better agreement with the data: the texture Bj H performs relatively 
well for 023 while many points obtained from C IH fall in the 3er range of $13, and in every case very few realizations 
give the correct r (much less then in the normal hierarchy case). 

When we include the requirement for the three mixing angles and r to be simultaneously into their 3cr confidence 
levels, the fractions of successful trials is almost the same: 



{B[ ,B^ H , B z , B{ H ,C 1H } = {0.0008, 0.0007, 0.0008, 0.0008, 0.001] 



(16) 
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B^ H (dashed line) and C IH (thick dashed line). The vertical black lines indicate the 3a allowed region on r and sin 2 #13. 
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As for the NH, we investigated the correlations among 5,i and 02,3, \an\ for the textures Bi and C (given the large 
2 rain , we do not need to study the Aj's). The qualitative behavior of the correlations for the B\ H textures follows 
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FIG. 8: Correlations among Si and £2,3 , |54 1 for the textures studied in the paper, for the case of IH. From above to below: 
Ai, A2, Bi, B2, B3, B4, C. 



that of the NH case since, as explained in SectHH the Si's describe equally well the NH and IH cases depending on the 
octant of 823- As expected from those same considerations, the parameters of the C IH are now less correlated than 
in the NH case, since the C mass matrix is more suitable to describe the inverted ordering of the mass eigenstates. 



V. CONCLUSIONS 

In this paper we have studied the compatibility of certain type of two-zero Majorana neutrino mass textures with 
the recent data in the neutrino sector. Differently from what has been extensively studied in the literature, we used the 
neutrino data to get information directly on the parameters entering the mass textures, with the aim of understanding 
how naturally a given texture is able to reproduce the value of the mixing angles and mass differences. To achieve 
this aim, we have used two different estimators: the minimum of the x 2 , built to extract the values of the matrix 
elements which better reproduce the experimental data, and the variable cIft, defined in eq.©, which is essentially 
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a measure of the fine-tuning among the parameters. 

We performed such an analysis for both normal and inverted ordering of the neutrino masses. In the first case, 
we observed that the patterns Bi and Ai fit the data with a very good x 2 > although the textures A 12 seem to be 
more favored. In fact patterns Ai also show a relatively small dpr and a less pronounced hierarchy among the matrix 
elements, features not shared by the patterns Bi. Texture C gives the largest (1ft and in addiction it needs #23 very 
close to the maximal value, that currently seems disfavored. Summarizing, textures Bi and C can be classified as 
less natural with respect to Ai because of the larger fine-tuning parameter dpT and hierarchies among the matrix 
elements. 

Taking the central values of the obtained dj's, we have derived approximate expressions of the textures in terms of 
the Cabibbo angle; extracting random 0(1) coefficients, we have seen that such textures show a similar agreement 
with the data as for the original patterns. Again, the textures derived from A t give a larger fraction of predictions for 
the mixing angles and the ratio r compatible with the data at the 3cr level; the correlation among the matrix elements 
is not as pronounced as in the other cases (especially for C, where the agreement to the experimental data is only 
possible at the prize of a huge conspiracy among the a, parameters). 

In the inverted case, the situation is more puzzling; in fact, from the fit procedure we see that the pattern C is less 
affected by fine-tuning and hierarchy problems than the Bi textures (the A\^ patterns are in much worse agreement 
with the data and have not been further studied); the approximate texture obtained from C produces large mixing 
angles (with its matrix elements not suffering from strong correlations) and can be considered particularly suitable to 
describe large 8i 3 . 
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